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ABSTRACT 

While the non-thermal radio through at least near-infrared emission in the hard state in X-ray 
binaries (XRBs) is known to originate in jets, the source of the non-thermal X-ray component is 
still uncertain. We introduce a new model for this emission, which takes into account the transient 
nature of outflows, and show that it can explain the observed properties of the X-ray spectrum. Rapid 
radiative cooling of the electrons naturally accounts for the break often seen below around 10 keV, and 
for the canonical spectral slope F v oc v' 1 / 2 observed below the break. We derive the constraints set by 
the data for both synchrotron- and Compton-dominated models. We show that for the synchrotron- 
dominated case, the jet should be launched at radii comparable to the inner radius of the disk (~few 
100 r s for the 2000 outburst of XTE J1118+480), with typical magnetic field B > 10 6 G. We discuss 
the consequences of our results on the possible connection between the inflow and outflow in the hard 
state of XRBs. 

Subject headings: radiation mechanism: non-thermal — stars:winds, outflows — X-rays: binaries — 
X-rays: bursts — X-rays: individual: XTE J1118+480 



1. INTRODUCTION 

The flat/inverted radio spectrum seen in the hard 
state of black hole X-ray binarie s (XRBs, e.g., 
iHvnes et qZ~ll2000t iFender ll200ll 120061 ). is naturally ac- 
count ed for by self-absorbed synchrotron emission from 
iets (iBlandford fc Konigl I Il979t iHiellming fc Johnston 1 
Il988t Falckc fc Biermann I I1999D , analogous to the com- 
pact radio cores of active galactic nuclei (AGN). Al- 
though generally too small to resolve, th e compact jets 
have been imaged in the case of Cyg X-l (jStirling et al, I 
.2001]), confirming the origin of the emission. In con- 
trast, the origin of the hard, non-thermal X-ray com- 
ponent is more uncertain. On the one hand, it has 
long been known that accretion-based models (inverse 
Compton emission from the inner parts of the disk) 
can have a si gnificant contribution to the hard X-ray 
ectr um (e.g."lTitarchuk ifliM M; Magdziarz & Zdzi arski I 



iGierlinski et al. I Il997t lEsin et al 
Poutanen I 119981 : 
2007i r On the 
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iCadolle Bel et al. I 120061 : lYuan et aU 
other hand, it is also possible that 
jets can contribute signifi cantly to the observational 
properties at these bands dMirabel fc Rodriguez I Il994l: 
Markoff et aT\ I2001L 120031 120051 : iBosch-Ramon et al. 
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2006; Parcdcs et al 120 06: Kaiser 200 9: iGupta et al. 
200a Kvlafis et al. 1 120081: iMaitra et al. 1 12009ft . Indeed 



the emission from jets in other objects, such as AGN, 
can span decades in spectral frequency, from the radio 
all the way to the X- and 7-ray bands. 

While both scenarios can pr ovide equally good fi ts to 
the X-ray spectrum (see, e.g., iNowak et al. ll20lTl ). the 
thermal corona inverse Compton (TCIC) classes of mod- 
els do not generally connect well to the need for strong 
magnetic fields and outflows in the same region, implied 
by the necessity of launching collimated jets. Because 
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of the very high temperatures in the inner disk regions 
(~ 1 keV: iShakura fe Sunvaevl 119731 ) . the hard X-ray 
emitting plasma will be mostly ionized, and thus be sub- 
ject to electromagnetic forces and in potentially near rel- 
ativistic motion. This motion together with the size of 
the dominant emitting region implies a characteristic dy- 
namical time that can play an important role in the ra- 
diative properties of the system, and leave its mark in 
the spectrum. 

Additional clues come from the spectral behavior at 
the X-ray band. In several outbursts, the X-ray spectrum 
shows a break between a few to ~ tens of keV, where the 
spectral index steepens from a ~ 0.5 (where F v cx v~ a ) 
to a ~ 10. Such a break can be seen in numerous objects. 
For example, it is seen at 2 keV in the 2000 outburst o f 
XTE J1118+480 (|Hvnes et aZ. 1120001: lEshTet al. 1120011) , 
at < 40 keV in the 20 00 outburst of XTE J1550-564 
([Rodriguez e t al. 1 12003ft. at ~ 70 k eV in the 2002 out- 
burst of GX 339-4 (|Homan et al. I [20051) . at 40 keV 



in th e 2005 outburst of GRO J1655-40 (|Joinet et al. I 
2008) , at ~ 40keV in the 1991 outburst of Cyg X-l 
(|Gierlinski et al. 1119971 ) and at ~ 7 keV in the 2002 out- 
burst of GRB 1758-258 ([Soria et al. |[20Tll) . to name only 
few examples. These breaks are most pronounced in the 
'hardest' state. Often, these breaks disappear and the 
spectrum becomes softer with the increase of the lumi- 
nosity, and the tr ansition to the high/soft state (e.g., 
iDunn et al. ll20Tl . 

The existence of this class of spectral breaks is not eas- 
ily explained in the framework of current models. Several 
ideas include Compton reflection by cold matter, which 
is ex pected to harden the X-ray spectrum above ~ 10 
keV dLightmn fc White1l98alHaardt fc Maraschi 1119931 : 
Magdziarz & Zdzia rskilll995l ). Anoth er idea is metal ab - 
sorption by partially ionized gas (e.g. lEsin et iDEooil). 
Alternatively, the break may result from combination of 
synchrotron spectrum at low energies, and synchrotron 

3 Often the data at these energies are noisy, hence the errors in 
the values of a after the break are fairly large. 
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self-Compton at higher energies (e.g., Markoff et al. I 
[200llHoman et al. Il2005h . 

While the above mentioned suggestions may indeed 
play a role in producing the break, we note that the low 
energy spectral index, a ~ 0.5, may hold the key to 
understanding the emission. Indeed, this spectral index 
is a natural outcome of both synchrotron and inverse- 
Compton emission processes, provided that the radiat- 
ing particles cool on a time scale much shorter than the 
dynamical time scale. This fact together with the need 
to account for plasma bulk motion thus motivates us to 
suggest an alternative model for the hard X-ray spectrum 
in hard state XRBs. 

As we show below, the observed X-ray spectral prop- 
erties of many of these objects can follow directly from 
the acceleration of electrons at some radius r. While 
it is clear that particles are accelerated to high energies 
in many astronomical objects, the details of this pro- 
cess are still highly uncertain. From a theoretical per- 
spective, it has long been argued that Fermi accelera- 
tion in shock waves (lBelT1[l978t IBlandford fe Ostriker I 
Il97l IBlandford fc Eichler I I1987h. or possib ly in mag- 
netic reconnection layers (e.g.. fKirkTI 2004 and refer- 
ences therein) is a plausible mechanism. Many works 
study this phenomena in b oth the Newtonian regime 
(e.g., lAmato fc Blasi 1 120051). an d the ultra-relativistic 
regime (e.g.. ISpitkovskvl 120081: [S ironi fe Spitkovskvl 
120111 ). These works suggest that indeed electrons are 
accelerated in shock waves into a power law distribution 
above a characteristic Lorentz factor, "f m . However, it 
is found that in both regimes only about ~ 1% of the 
particles are accelerated to a power law, the rest main- 
tain a low energy Maxwellian distribution. A similar 
result was found in studies of particle acceleration in su- 
pernovae remnants (SNRs; see lEllison et al. 1 12007). The 
emerging theoretical picture is therefore that particles 
are heated as they cross the shock wave, and their spec- 
tra is composed of a relativistic Maxwellian distribution 
with a small fraction accelerated to a power law distr i- 
bution at higher energies (see also iCaprioli et aTl^OlOj ). 

On the other hand, phenomenological works of spectral 
fitting in various astronomical ob jects (e.g., astrophysi - 
cal jets, supernovae remnants - ILazendic et al. ||2004D 



and gamma-ray burst (GRB) a fterglows (|Galama et al. I 
Il998t iWiiers fe Galamal Il999f ) suggest that the power 
law tail of the electron distribution dominates. Thus, 
there is a discrepancy between the phenomenological 
works and the theoretical predictions. 

Regardless of the exact details, following the acceler- 
ation the electrons radiate their energy via synchrotron 
emission and inverse-Compton scattering of both the disk 
(thermal) UV/X-ray photons and locally produced syn- 
chrotron photons (synchrotron self-Compton; SSC). As 
the electrons radiate their energy they cool. I ntere stingly 
enough, it was shown by iPe'er fc Casella I ([2009) , that 
the only requirement needed in order to obtain the uni- 
versal spectrum observed below the break, F v oc j/ -1 / 2 , is 
that the cooling time scale is shorter than the dynamical 
time scale (the particles are in the "fast cooling" regime) . 
Thus, this result is independent on the uncertain details 
of the acceleration process. 

The spectral shape above the break does depend 
on the details of the acceleration process. Due 



to the uncertainty mentioned above, in the calcula- 
tions presented below we consider both scenarios, the 
theoretically-motivated scenario in which most accel- 
erated particles have a Maxwellian distribution, and 
the phenomenologically-motivated scenario, in which the 
spectral shape of most accelerated particles is a power 
law. 

Fast cooling of a power-law distribution of electrons 
leads to a spectral index F v oc v~ p l 2 above the spectral 
break, where p is the power law index of the accelerated 
electrons. Thus, if most electrons are accelerated to a 
power law distribution, then the observed spectral slopes 
at high energies, a ~ 1, are therefore consistent with the 
nearly universal power law index p — 2.2 ±0.2 observed 
in many astronomical objects. As we claim here, the 
derived conditions at the emission site, such as emission 
radius and magnetic field strength are most naturally 
accounted for in turbulent outflows, as would be expected 
for a region of jet launching. 

If, on the other hand, most accelerated electrons have 
a Maxwellian distribution, than a power law spectrum 
above the break is not expected. However, as we show 
below, the finite width of the Maxwellian distribution 
makes it possible to fit the observed spectra in this sce- 
nario as well. In this case, the existing data does not 
provide information on the power law index of the most 
energetic electrons. 

We therefore present here a new jet model that can 
account for the X-ray properties of many X-ray bina- 
ries. By solving the kinetic equations that govern parti- 
cle and photon populations self-consistently, we are able 
to reproduce the main observational properties of the 
X-ray spectra. In particular, we suggest a natural ex- 
planation for the observed a — 0.5 to a ~ 1 spectral 
break as resulting from rapid radiative cooling of the en- 
ergetic particles. Our num erical calculations are ba sed 
on the model developed by IPe'er fc Waxmanl (j2005l ) in 
the study of the prompt emission of GRBs, that was fur- 
ther mod ified in the study of the radio properties of X-ray 
binaries (|Pe'er fc Casella 1I2009D . Here, because we focus 
on the X-ray properties of XRBs, various physical pro- 
cesses that were not required for radio emission studies 
are now included, such as Compton scattering and pair 
production and annihilation. 

This paper is organized as follows: In [21 we present the 
basic properties of our model. In particular, we focus on 
the rapid cooling that leads to the universal F v oc v~ x l 2 
spectra below the spectral break. We compare our model 
to the data in Sj3j before summarizing and discussing our 
results in view of the recent observations in [j4] 

2. MODEL ASSUMPTIONS AND BASIC PROPERTIES 



Sim ilar to earlier j et models (e.g.. lMarkoff et aI~ll200lL 
2005; iMaitra et al\ I2009T) . we assume that the initial 
gas which is transferred from a mass-losing secondary 
forms a geometrically thin , optically thick cool disk 
(jShakura fc Sunyaev I [1973). Emission from the disk 
dominates the emission at the extreme-ultraviolet and 
soft X-ray bands. Thus, fitting the data at these energy 
bands provides an estimate of the thermal luminosity L t h 
and the temperature Td of particles and photons in the 
inner part of the disk, and is often seen as a proxy to 
the bolometric luminosity of the system. The inner disk 
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radius is then given by rj„ ~ (Ltft,/27ro's) 1 / 2 TJ~ 2 , where 
as is the Stefan-Boltzmann constant. 

There is still some debate over whether the data from 
hard state XRBs implies a truncated inner disk radius or 
not, with often conflicting results from continuum fitting 
of the thermal accretion disk, and the modeling of the 
potentially broadened iron fluorescence line at 6.4 keV. 
Several groups claim an inner radius of typically Ti n ~ 
10 2 - 10 r s , where r s = 2GM B h/c 2 is the Schwarzschild 
radius of a non- rotating black hole (lEsm et al. 1 12001 
iTomsick et 0171120091 IDone fc Diaz Trigo ll2010ir ~On the 
other hand, several recent works argue t hat the data im- 
ply a lower inner radius, ~ 10 — 10 15 r s , ((Reynolds et al. I 
12010ft . or even agai nst the need for a recessed disk to ex- 
plain the data (e.g. JReis et al. ll2009HRevnolds fc Miller I 
12010ft . 

Regardless of the exact radius where it occurs, gen- 
erally we do expect a transition of properties in the 
inflow/outflow of the accreting plasma. Once the gas 
approaches the disk inner radius r^, it is expected to 
become geometrically thicker and optically thin, and 
also radiatively inefficient. Several variations of the 
original advection-dominated accretion flow (ADAF; see 
iNaravan fc McClintockl I2008L for a recent review) sce- 
nario, more generally referred to as radiatively ineffi- 
cient accretion flows (RIAFs) now exist. Because RIAFs 
are likely necessary to support jet launching (see, e.g., 
ILivio et al. Ill999l ; (Meier ll200l[) .' we here assume such a 
transition is likely. However, since emission from RIAFs 
has been ex tensively studied in the past (e.g.. lEsin et al. I 
119971 [200ll ) . we currently focus on the jets. While we as- 
sume here that the jets dominate the emission^] , we note 
that a possible contribution from the RIAF may alter 
our final results. 

We thus assume that at some radius rj, part of the 
plasma is ejected outward from symmetric nozzles into a 
jet outflow. In the calculations presented, rj is taken as 
an independent variable, that is not necessarily equal to 
ri n . During the ejection process, the electrons are accel- 
erated. Such an acceler ation can result from Fermi accel- 
eration in shock waves (Ellis on et a7~lll990t iSpitkovskv I 
2008) , or perhaps in reconnection layers of stron g magne- 
tized outflows (e.g., Drcnkh ahn fc Spruit II2002I ). We pa- 
rameterize the current uncertainties about the details of 
the acceleration process by simply assuming that a frac- 
tion £ p i of the electrons are accelerated to a power law 
distribution with power law index p above a characteris- 
tic Lorentz factor 7,„, likely set by the thermal tempera- 
ture of the innermost accretion flow. The remaining elec- 
trons (a fraction 1 — £ p ;) assume a relativistic Maxwellian 
distribution, with temperature 6 = kT/m e c 2 — 7 m /2, 
ensuring smooth connection with the power law at higher 
energies. 

In order to connect the acceleration microphysics to the 
system energetics, we further assume that a fraction e e of 
the jet kinetic luminosity is dissipated away, carried 
by the energetic electrons. This assumption therefore 



Similar assumptions were taken by, e.g., lYuan et a£\ ([200l k 
lYuan fc Cml (2003) 



results in the constraint (e.g.. IPe'er fe Waxman I [2004T ) : 

'-j::|i£ (=l ;:" - 

Here, e max (e m i n ) is the maximum (minimum) energy of 
the accelerated electrons . The maximum energy is cal- 
culated by equating the acceleration time to the cooling 
time (e.g., via synchrotron emission), and the minimum 
energy, for power law index p = 2 is calculated iteratively, 
as e m j n = 7 m m e c 2 . Following the calculation of 7 m , it is 
corrected iteratively to account for the fraction 1 — £ p i of 
the electrons that assume a Maxwellian distribution, so 
that the energy density in the combined (Maxwellian + 
power law) electron population exactly equals a fraction 
e e of the energy density in the jet. 

The particles propagate at bulk velocity /3jC inside the 
jet, during a (comoving) dynamical time td yn — lj r jl PjC, 
where jj = (1 — /3 2 )~ 1//2 is the Lorentz factor associated 
with the bulk motion of the plasma. During this time, a 
continuous injection (and acceleration) of particles at the 
base of the jet is assumed. We assume the existence of a 
steady magnetic field of strength B inside the jet. While 
in some works the magnetic energy density is described 
as a fraction eg of the jet energy (analogue to the defi- 
nition of e e ), we note that as the source of the magnetic 
field may be attached to the inner parts of the disk, it is 
possible that es > 1. Therefore, here we consider B to 
be a free parameter. As a result, once introduced into the 
jet, the particles lose their energy via radiatively produc- 
ing energetic photons. This energy loss results from syn- 
chrotron emission and by inverse Compton scattering of 
both the thermal disk photons as well as the synchrotron 
emitted photons (SSC). 

The existence of a characteristic physical scale rj (cor- 
responding to characteristic time td yn ) implies that the 
physical quantities (e.g., the magnetic field B or the en- 
ergy density in the photon field) do not vary significantly 
over this scale. In particular, we neglect adiabatic energy 
losses over this scale, as its contribution to the electrons 
cooling is (by definition) much w eaker than the radiativ e 
cooling (see further discussion in IPe'er fc Caseflal l2009') . 

The radiative cooling time of electrons at energy E = 
7m e c 2 is given by (in the Thompson regime): 

, E jm e c 2 

t C ooi- p {mcaTl 2 {uth + UB y W 

where u th = L th /nr 2 f3jC^j and u B = B 2 /8n are the 
(comoving) energy densities of the thermal photons (as- 
suming cylindrical geometry) and the magnetic field, re- 
spectively. We denote by j c the Lorentz factor of the 
electrons for which the cooling time is comparable to the 
dynamical time, 7 C = (m e c 2 ) / (A/3)ca T (u th + u B )t dyn . 
For 7 C <C 7 m , the electron energy distribution in the 
range j c <C 7 <C 7 m reaches a steady state, which is 
calculated by solving the rate equation, dn e i(^,t)/dt — 
((5/(97) [ n ei{l)d"f / dt] = 0. Since the power emitted by 
both synchrotron radiation and Compton scattering is 
Psynjc dj/dt oc 7 2 , the steady state electron dis- 
tribution at the range 7 C < 7 C j m is ?i e ;(7) oc 7~ 2 . 
Above 7 m , the steady particle distribution (assuming an 
initial power law) is n e /(7) oc 7 _<J)_1 ). Therefore, the ob- 
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tained spectral shape from either synchrotron- or inverse 
Compton-dominated cooling is 

where frv m — £ m i n is the typical energy of photons emit- 
ted by electrons at 7 m . 

Equation |3] is valid if most electrons are accelerated to 
a power law. If -C 1, then above the break at v m a 
smooth variation in the spectrum is expected, resulting 
from the finite width of the Maxwellian distribution (see 
below) . 

The solution obtained in equation [3] is valid only as 
long as 7 C — 1 > 1 (and 7 C <C 7™)- If the cooling is ex- 
tremely rapid, the electrons become sub-relativistic on a 
time scale much shorter than the dynamical time. Un- 
der these conditions, the sub-relativistic (cold) electrons 
interact with the low energy (disk) photons via multi- 
ple Compton scatterings, altering their original distribu- 
tion (for discussion o n this effect in the contex t of GRB 
prompt emission, see lPe'er et al. I [20051 [2006'). Thus, a 
sharp cutoff of the disk photons can be used to constrain 
the optical depth, hence discriminate between jet and 
accretion-based models. See further discussion in £|4] be- 
low. 

3. CONSTRAINTS SET BY OBSERVATIONS 

The broad band observations of the XRB hard state 
contain a wealth of data that can be used to provide 
strong constraints on the free model parameters. The 
rapid cooling of the electrons implies that nearly 100% of 
the energy given to the electrons is converted to emission, 
mainly at the X-ray band. Thus, a constraint on the jet 
kinetic luminosity as well as the fraction of energy used 
to accelerate the electrons is Lnt — Lke e , where Lnt is 
the luminosity in energetic photons (with energies above 
that of the disk photons). The break energy at few keV 
corresponds to £ m i n , and, for £ p j < 1, the spectral slope 
above the break provides a direct indication of the power 
law index p of the accelerated particles via eq. El 

A large optical depth to scattering by the electrons 
and the produced pairs will result in a modification in 
the emission from the inner parts of the disk. Thus, the 
sharp cutoff in the disk emission implies that the optical 
depth to scattering cannot be much larger than unity, 
r 7e P$ 1- An order of magnitude estimate of the optical 
depth is obtained by assuming cylindrical jet geometry 
with scale rj. The comoving number density of injected 
electrons is n' el = Lkl^r 2 m p c 2 ^j 2 , resulting in optical 

depth T 7e ~ ~ijTjn' e icrT = Lnt&t /Trrje e m p c 3 jjf3. 

The optical depth is modified by the existence of pairs. 
An estimate of the optical depth in the vicinity of pairs 
can be done as follows (a detailed description will be 
presented in another work). For short dynamical times, 
pair annihilation is insignificant, and the number den- 
sity of pairs is n± ~ n 7 ( e>meC 2)n 7 ( E<meC 2)CCTT£<ij/n/4. 
For a flat spectral energy slope (yF v oc v°) over 
the energy range £ m in _ £max (as is expected for elec- 
trons power law index p not much different than 2), 
the number density of photons at the different energy 
bands are n 7 ( e>mcC 2) ~ u NT /m e c 2 log(e max /e m in), and 
™ 7 (e<m c c 2 ) - ™ 7 ( e >m c c 2 ) x (jn e c 2 /e min ) , where u NT = 
Lnt I TTr 2 c/3^f 2 is the energy density in the non-thermal 



component. Combined together, the optical depth for 
scattering by the pairs is 



t- , ~ ( L nt ctt \ a. m e c 2 

T> e ~ \2^Fft^fn^ ) ^3 £ mi »log(£ max /e min ) 2 

~ i ( lnt \ 2 ( rj v 2 

- Jff* \5xlO- 3 L E J \Wr s J 



In estimating the optical depth in equation |4j we used 
parameters that charact erize XTE J 1118+4 80, such as 
distance d = 1.72 kpc (jGelino et al. Il2006h . We esti- 
mate Mbh — 7M(r0, resulting in Eddington luminosity 
L E = 8.75 x 10 38 ergs _1 . Since log(£ max /e min ) w 10, for 
e min ~ fewkeV, rn e c 2 /£ min log(£ max /£ min ) w few. Thus, 
the requirement r 7 e ± < 1 implies, for mildly relativistic 
jets (say, j3j ~ 0.4), rj > 10 15 r s . Note that the charac- 
teristic jet scale 77 and the jet (non-thermal) luminosity 
Lnt are considered as independent variables, measured 
in units of r s and L e ■ 

The break seen at few keV is attributed to emission 
from electrons at 7 m , either via synchrotron emission or 
Comptonization of the disk photons. As these two sce- 
narios lead to different constraints on the model param- 
eters, we treat each one separately. 

In the calculations below, we scale the values of the 
model parameters to the reference values relevant for the 
2000 outburst of XTE J1118+480. This is done for ease 
of comparison with the numerical results presented in 
33.31 As the discussion is general and the dependence on 
the spectral breaks and luminosities are given, general- 
ization to any source is readily obtained. 



3.1. Synchrotron- dominated model 

Electrons with Lorentz factor 7 m radiate synchrotron 
photons at characteristic observed energy = 
(3/2)h(qB/m, e c)^V = 1.75 x 10~ 8 B^V eV. Here, q 
is the electron's charge, and V = [7^(1 — j3j cosO)]^ 1 > 1 
is the Doppler shift (8 is the angle to the line of sight). 
The requirement that £ m j n ~ 2keV thus constraint the 
values of i^^J!! A further constraint on B is added by 
the physical requirement that e e < 1, which provides an 
upper limit on the value of j m via equation [TJ 

The second requirement is that the cooling break, v c , 
which is the characteristic frequency of synchrotron emis- 
sion from electrons having Lorentz factor 7 C is obscured 
by the disk photons, hence e° b = hvf> < 0.3 keV. Fi- 
nally, we require that synchrotron emission is the domi- 
nant cooling mechanism, which translates into m_b S> Uth- 
Combined with the requirement on the optical depth (eq. 



5 While IMcUlintock et al. I 12001]) estimated M BH = 6.1M , 
IGelino et al. I ^2009 ) obtained Mbh = 8.5M . 

6 For mildly relativistic jets, 7^ > 1, the Doppler shift is not 
strongly constrained. 
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0}, one obtains the set of constraints 

1. a. e°£ ~ 2keV : B 7l 2 n ~ 1.2 x lO 11 !?" 1 G 
1.6. e e < 1 : B > 3.4 x 10 4 2? _1 G 

2. e° b <0.3keV: 

B > 4.3 x 10 4 (ft-/7 i ) 2/3 (r- i /10r s )- 2 / 3 X' 1 /3G 

3. u B » itt/, : 

B»1.3x UPpT 1 ^- 1 (Lth/3.2 x l(r 3 £ £ ) 1/2 x 
(r 3 -/10r a ) G 
4- r 7 . e± < 1 : fo/lOr.) > ipf'^lj 1 ' 2 

In equation [5j the thermal (disk) luminosity is normal- 
ized to 3.2 x 1CT 3 L B , as is in the 2000 outburst of XTE 
J1118+480; see further discussion in £13.31 below. 

One thus obtains strong constraints on the value of 
B > 10 6 G and the characteristic scale rj > 10 15 r s . 
While the value of e e is not strictly limited from below, a 
high value of e e < 1 is generally preferred, since lower val- 
ues of e e implies, via eq. [5] (l.a.) stronger magnetic field, 
which may be challenging to produce. High value of e e 
implies, in turn, high radiative efficiency - namely, that 
the jet kinetic luminosity is similar to the non-thermal 
luminosity seen. We further note that there is no con- 
straint on the value of £ p z set by the data below e^. 

3.2. Compton- dominated model 

Alternatively, the break seen at ~ few keV may be at- 
tributed to inverse Comptonization of the disk photons 
by the jet electrons. As the angle-averaged Comptonized 
energy is e m — (4/3) 7 ^£i n where ~ 0.1 keV is the en- 
ergy of the disk photons (using the value from the 2000 
outburst of XTE J1118+480), in this scenario the charac- 
teristic electrons Lorentz factor is 7 m < 10. This result, 
in turn, implies e c >3x 10" 3 , and L k < Le- 

In this scenario, the cooling is dominated by Comp- 
tonization of the disk photons, u t h S> ws. Finally, one 
requires the cooling time to be shorter than the dynami- 
cal time, which implies j c > 1. Combined together, these 
conditions give 

1. ~ 2keV : e e > 3 x 10" 3 

2. 7c - 1 : (r/10r.) ~ Wflj 1 (W3.2 x IQ- 3 L E ) 

3. u B < u t h ■ 

S«1.3x laPpr 1 '*^ 1 (Lth/3.2 x lO- 3 ^) 172 x 
(rj/10r s ) G 
4- r 7 , e± < 1 : fo/lQr.) > l/3J 3/2 ^ 1/2 . 

(6) 

We stress again that while the derived constraints on 
the parameters values are obtained using the data from 
the 2000 outburst of XTE J1118+480, the numerical val- 
ues for other sources are easily obtained by inserting the 
values of the break energies and the thermal and non- 
thermal luminosities in equations [SJ [6l 

3.3. Numerical results 

The analytical constraints derived in equations [5j [6] 
provide good approximations for the required condi- 
tions from jet model needed in fitting an observed spec- 
trum. However, due to the complexity of the prob- 
lem, in particular the possible creation and annihila- 
tion of a large number of pairs, a numerical approach 



is needed in order to validate the analytical approxi- 
mations and provide more accurate constraints. Here, 
we use a ti me-dependent mod e l, bas ed on the code de- 
veloped by IPe'er fc Waxman I ((20055. This code solves 
self-consistently the kinetic equations that govern the 
time-dependent energy distribution and radiation from 
particles inside the jets, following all relevant radiative 
processes. These include synchrotron emission, Compton 
scattering, pair production and annihilation that follow 
particle acceleration at the base of the jet. An example of 
the numerical results are presented in Figures [TJ [5] on top 
of the data from the 2000 outburst of XTE J1118+480. 
We stress that no attempt was made to obtain the statis- 
tical best fit to the data, but only to demonstrate the pos- 
sibility of reproducing the observed X-ray spectral slopes 
in the two model scenarios considered. 

In Figures [T] and [H we show the results of our explo- 
ration of two scenarios, strong particle acceleration (£ p ; = 
1) and weak particle acceleration (£ P 2 = 0.01), respec- 
tively. In both figures, the solid (blue) line represents the 
numerical result for the synchrotron-dominated model. 
We chose representative parameters consistent with a 
mildl y-relativistic jet (/?,- = 0.4), s imilar to the value cho - 
sen in lFalcke k, Biermann I ()1995l ); IMarkoff et al. I(|2001D . 
The disk is fitted with a multi-color black body compo- 
nent, with L t h — 0.0032Ls, and inner disk tempera- 
ture Ti n = 5 x 10~ n erg, giving an inner disk radius 
of ri n /r s — 326. We chose parameters values consistent 
with the constraints set in equation [SJ such as magnetic 
field B = 3 x 10 6 G and rj /r s = 200. The value of e e is 
close to equipartition: e e = 1.0 in Figure [T] and e e = 0.5 
in Figure [21 The resulting spectrum is very similar to the 
X-ray spectrum of XTE J1118+480, for L k = 0.0045L B 
(Figure [J) and L k = 0.005i_e (Figure E} . In the results 
presented in both figures, the power law chosen above 
7m is p = 2.0, consistent with Fermi acceleration mod- 
els. We note that in the £ p i = 0.01 scenario presented in 
Figure [21 the high energy part of the spectrum is fitted 
with particles at the peak of the Maxwellian. Reasonable 
fits are obtained due to the finite width of the Maxwellian 
distribution. In this scenario, however, the exact value 
of p cannot be determined from the data. 

For these parameter values, the flow is (by definition) 
magnetic dominated: the energy density in the magnetic 
field, ub = B 2 /8ir ~ 3 x 10 11 erg cm -3 is much larger 
than the kinetic energy density, u k — Lb/nr 2 /3jCjj ~ 
5 x 10 8 erg cm" 3 . 

The result of a Compton-dominated model is presented 
by the dashed (green) line in Figures Q] and O The values 
of the disk parameters and the jet velocity are similar to 
the ones used in the synchrotron model. The values of the 
free parameters of the jet model chosen, rj/r s — 25 and 
e e = 0.029 (Figured]), e e = 0.01 (Figure [2]) reproduce 
the data well for L k = 0.15L E (Figure [Q), L k = 0.2 
(Figure [2]) , and weak magnetic field. Here we chose 
B = 10 3 G, which implies matter dominated outflow 
- UB/u k ~ 3 x 10~ 8 . As explained above, the exact 
value of the magnetic field is unimportant as long as the 
constraints in equation [5] are fulfilled. In this scenario as 
well, for ^ p i — 1 a power law index p = 2.0 results in very 
good fit to the data. 

4. SUMMARY AND DISCUSSION 
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Fig. 1. — Example of the numerical results, on top of the data of 
the 2000 outburst of XTE J1118+480. The solid (blue) curve rep- 
resent the synchrotron-dominated model, while the dashed (green) 
curve represents the inverse Compton-dominated model. A power 
law- dominated distribution of accelerated electrons (£ p ; = 1) is 
assumed. The values of the other model parameters are given in 
the text. No attempt was made to statistically fit the data, but 
rather this plot demonstrates the ability of this transient jet model 
to reproduce the basic features of the observed X-ray spectra. The 
dotted (red) line is (multi-color) black body component, originat- 
ing from the accretion disk. These photons (together with the syn- 
chrotron photons) serve as seed photons for Compton scattering 
by the energetic electrons. 




Obs. energy [eV] 

Fig. 2. — Similar to Figure[T] but with only a small fraction (£ p ; = 
0.01) of particles accelerated into a power law, the rest maintain a 
Maxwellian distribution. As explained in the test, reasonable fits 
are obtained in this scenario as well, due to the fact that in the 
fast cooling regime, a universal spectrum F v oc v 1 / 2 is produced 
regardless of the initial electrons spectra, and that at high energies, 
the finite width of the Maxwellian enables fitting the data. In this 
scenario, the fit to the data does not uniquely determine the value 
of power-law index p. 

In this paper, we present a new jet-dominated model 
that is able to reproduce the main spectral properties 
seen at the X-rays of many XRBs in the hard state. Our 
key motivation is the spectral break which is often seen 
at ~ few keV - few tens keV , and (within these sources) 
the nearly universal spectral slope, F v oc v 1 / 2 observed 
below this break. This spectral slope is a natural out- 
come of emission from electrons whose energy distribu- 
tion is determined by rapid radiative cooling, following 
acceleration at the jet base. The rapid cooling can result 
from either synchrotron emission or Comptonization of 
the disk photons, or as likely, a contribution from both 
processes. We derive in equation [5j [6] the required con- 



straints on the free model parameters for both these sce- 
narios, and demonstrate the resulting spectrum in Fig- 
ures [T] and [2] 

The key difference between our model and earlier jet 
models is that here we self-consistently consider the tem- 
poral variation of the particle distribution due to the 
radiative cooling, and consider those variations in the 
spectral calculations. Thus, we are able to make a clear 
separation between the acceleration process and the cool- 
ing processes; namely, we do not assume that the en- 
ergy lost by the electrons is necessarily fully replenished 
by any heating source. Indeed, in recent years it has 
become clear that this decoupling between acceleration 
and cooling may hold the key to understanding the spec- 
tral properties of XRBs. As a result, in recent years 
several numerical models which consider this separation 
have been constructed in the study of XRB s spectra (e.g., 
IBelmont et~aT1 12008t [Maitra et aH 12009ft . In essence, 
these models are very similar to the numerical model 
used here. A key advantage of our model is its ability of 
solving the kinetic equations over many orders of magni- 
tude in time and energy scales, which makes it ideal for 
producing broadband spectra. 

As indicated in Figures Q] and [51 both synchrotron- 
and Compton-dominated scenarios can in principle re- 
produce the X-ray spectrum of a typical source show- 
ing the predicted break. The values of the free model 
parameters are, however, significantly different in the 
two scenarios. While the Compton dominated sce- 
nario requires r/r s < 10 1,5 , ~ Le and a rela- 
tively weak magnetic field, the synchrotron dominated 
scenario better fits the data with r/r s > 10 2 , Lk < 
10~ 2 Le and a strong magnetic field, B > 10 6 G. In 
our opinion, at least for XTE Jl 118+480, the data fa- 
vor a synch rotron-dominat e d scenario, similar to the con - 
clusions of iMarkoff et al. I (|2001[ ); iMaitra et aT\ (|2009ft . 
First, the characteristic scaling rj is similar to the inner 
(truncated) disk radius, r;„. Thus, obviously, this radius 
marks a physical transition in the properties of the flow, 
which can result from development and launching of jets. 
Second, the non-thermal radio spectra observed in the 
hard states provide clear indication for jet-synchrotron 
emission. Emission at the radio band is expected once 
the plasma reaches a scale much larger than rj (and with 
corresponding much weaker magnetic field) . As the mag- 
netic field is expected to decay along the jet (due, e.g., to 
Poynting-flux conservation) , unless there is a generation 
of strong magnetic fields at some scale r 3> rj, a strong 
magnetic field at the jet base decays along the jet, and 
thus the same field may be the source of the radio emis- 
sion, further along the jet. While studying the full con- 
nection between the X-ray spectrum and the radio prop- 
erties is left for future work, we note that such a connec- 
tion may exist, and may lead to very interesting cons e- 
quences (|Pe'er fe Casella 1120091 ; ICasella fc Pe'er II2009D . 

The numerical results presented in Figures [1] and [2] 
show that the change in the spectral index (the spectral 
break) is gradual, extending over about an order of mag- 
nitude in energy. This result originates from the broad 
band nature of the emission processes. Thus, in partic- 
ular for scenarios in which the characteristic breaks 
and e° h are close to each other in energy, a gradual change 
of the slope is expected, providing a natural explanation 
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for the somewhat softer slope observed below the break 
in several objects. 

Similarly, a general trend seen in many objects is a 
softening of the spectrum with increase of the luminosity, 
marking the transition to the soft state (Migliari et al. I 
[20071 iDunn et al. ll20lTI ). This result can be understood 
in the framework of the model presented here as follows. 
As the luminosity increases, the conditions at the ac- 
celeration region are likely to vary. The strength of the 
magnetic field B as well as the acceleration efficiency (e e ) 
are likely decreasing, as is indicated by the lack of radio 
emission in the soft state. In both the synchrotron and 
the Compton dominated models, the break energy 
depends on e 2 e , and in the magnetized model also on B. 
Thus, as a result of the changing conditions, is shifted 
to lower energies, and the break is eventually obscured 
by the disk photons. At this stage, the X-ray spectrum 
represents the emission above the spectral break, which 
is softer. 

Another interesting result of this work is the ability 
to reproduce the observed spectra with acceleration pro- 
cesses that produce a power law index p ks 2 for the ac- 
celerated electrons, under the assumption that £, p i ~ 1. 
While a full theory of the acceleration process still does 
not exist, a plethora of observational evidence exists that 
indicates a canonical power law index of the accelerated 
electrons. These include evidence from AGN jets, su- 
pernovae remnants and GRB shock waves, all showing 
canonical index p > 2.0. Our results are thus consistent 
with these findings, and do not require a modification of 
this law. While the quality of the data above the spectral 
break in many objects is not very good, we note that of- 
ten a roughly flat power flux (vF v oc v ) is seen, further 
strengthening this result. 

While in recent years significant progress in un- 
derstanding particle acceleration in shock waves has 
been made, there is currently a discrepancy be- 
tween the theoretical expectation of j pi <C 1 (e.g., 
Amato fc Blasi 1120051 iSpitkovsky I I2008t fCaprioli etaU 
20101 IS ironi fc Spitkovskv 1 120111 ) and the phenomeno- 

logical fittings indicating highe r value of £„ ; ~ 1 

(e.g.. IWiiers fc Galama I ITqqI iFrail et al. \ l200fT: 
Freedman fc Waxman I 20011 : iPanaitescu fc Kumar I 



2002t iMarkoff et al. I I2005D . This discrepancy can 



perhaps be explained by the fact that the theory is not 
yet complete. Since this topic is outside the scope of 
this work, we have explored the consequences for both 
scenarios. 

Specifically, we have demonstrated that in the alter- 
nate scenario in which most accelerated electrons main- 
tain a Maxwellian distribution (£ p ; <C 1), good fits to the 
data can still be obtained. The high energy part of the 
spectrum is explained by emission from particles at the 
peak of the Maxwellian, although this scenario does not 
constrain the value of the power-law index p. 

While we show here that a jet model is able to re- 
produce the X-ray spectra, we note that possible con- 
tribution to the X-ray emission may arise from the in- 



ner parts of the disk (RIAF models). Observationally 
there are some new results that hint at an interplay of 
the two processes as a function of the total luminosity 
(I Yuan fc Cui I [20051 [Russell et al. I l20l3 IPlotkin et a,U 
120111: iGandhi et al. 1 120 111 ), where RIAFs may begin to 
dominate at the higher hard state luminosities. 

From a theoretical point of view, the conditions in the 
inner part of the inflow may be similar to the condi- 
tions in the inner pa rt of the outflow (see discussion in 
IMarkoff et al. Il2005l ). Hence, a clear separation between 
the inner disk contribution and the inner jet contribution 
may be difficult. Several methods may be used to sep- 
arate the disk and jet contributions. First, a clear con- 
nection between the X-ray emission and the radio/NIR 
emission, which must have a jet origin, could provide ev- 
idence for jet dominated model, provided that the phys- 
ical conditions at both emission sites (e.g., the magnetic 
field or the number of radiating particles) are matched. 
Such a connection already seems to be indicated by the 
observational papers cited above. 

A second method is based on the requirement that the 
optical depth to scattering of the disk photons needs to 
be not much larger than unity, otherwise the spectral 
shape of the disk will be altered. Thus, if the disk accre- 
tion time is significantly larger than ~ rj/c, this could 
enable multiple Compton scatterings between the elec- 
trons and the disk photons, which may introduce a differ- 
ent characteristic from inflow emission. A full treatment 
of this idea is left for a future work. Finally, the best 
way to discriminate between disk and jet models may 
be using X-ray polarization measurements. Both syn- 
chrotron emission and non-saturated Compton scatter- 
ing are highly polarized, hence jet emission is expected 
to be polarized. On the other hand, in RIAF models 
the polarization signal is expected to be weaker, due to 
averaging of the synchrotron signal from different parts 
of the disk, as well as contribution from free-free emis- 
sion from the outer accretion flow. We thus expect that 
proposed future X-r ay polarization miss ions, such as the 
POLARIX mission (jCosta et al. 1 120101 ) or the Gravity 
and Extreme M agnetism Small Explorer (GEMS) mission 
(Jahoda 2010) may play a significant role in determin- 
ing the relative contributions from the disks and the jets 
in producing the X-ray spectra. 
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